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O ABSTRACT 
O 
<N 

^ ■ We present HST/STIS long-slit spectroscopy of NGC 1976. Our goal is to 

Q ■ measure the intrinsic line ratio [0 III] 4364/5008 and thereby evaluate the elec- 

tron temperature (T e ) and the fractional mean-square T e variation (t\) across 
the nebula. We also measure the intrinsic line ratio [N II] 5756/6585 in order to 
estimate T e and t\ in the N + region. The interpretation of the [N II] data is not 
^ ■ as clear cut as the [O III] data because of a higher sensitivity to knowledge of the 

electron density as well as a possible contribution to the [N II] 5756 emission by 
(SI ■ recombination (and cascading). We present results from binning the data along 

the various slits into tiles that are 0.5" square (matching the slit width). The 
average [O III] temperature for our four HST/STIS slits varies from 7678 K to 
r£3 ■ 8358 K; t\ varies from 0.00682 to at most 0.0176. For our preferred solution, the 

average [N II] temperature for each of the four slits varies from 9133 K to 10232 K; 



Or 



t 2 A varies from 0.00584 to 0.0175. The measurements of T e reported here are an 
average along each line of sight. Therefore, despite finding remarkably low t A , 
we cannot rule out significantly larger temperature fluctuations along the line of 



sight. The result that the average [N II] T e exceeds the average [O III] T e confirms 
what has been previously found for Orion and what is expected on theoretical 
grounds. Observations of the proplyd P 159-350 indicate: large local extinction 
associated; ionization stratification consistent with external ionization by 9 l Ori 
C; and indirectly, evidence of high electron density. 
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Subject headings: ISM: abundances — ISM: atoms — ISM: H II regions — ISM: 
individual: NGC 1976 — ISM: individual: proplyd P159-350 

1. INTRODUCTION 

Most observational tests of the chemical evolution of the universe rest on emission line 
objects; these define the endpoints of stellar evolution and probe the current state of the 
interstellar medium (ISM). Gaseous nebulae (H II regions & Planetary Nebulae (PNs)) are 
laboratories for understanding physical processes in all emission-line sources, and probes for 
stellar, galactic, and primordial nucleosynthesis. 

There is a fundamental issue that continues to be problematic - the discrepancy between 
heavy element abundances inferred from emission lines that are collisionally excited (CELs) 
compared with those due to recombination/cascading, the so-called "recombination lines". 
Studies of PNs contrasting recombination and collisional abundances (Liu et al. 1995, Kwitter 
& Henry 1998) often find differences exceeding a factor of two. For NGC 7009, Liu et al. 
(1995) found that the recombination C, N, and O abundances are a factor of ~5 larger 
than the corresponding collisional abundances, which was found to be the case also for neon 
(Luo, Liu & Barlow 2001). For NGC 6153, Liu et al. (2000) found that C++/H+, N++/H+, 
++ /H + , and Ne ++ /H + ratios derived from optical recombination lines are all a factor of 
~10 higher than the corresponding values deduced from CELs. The discrepancy is even 
larger, a factor of ~20, for the Galactic bulge PN M 1-42 (Liu et al. 2001) and Hf 2-2, which 
has the most extreme abundance difference to date, a factor of 84 (Liu 2002). 

Studies of H II regions find similar behavior although the differences are generally smaller 
than is the case for PNs. For the Orion Nebula, Esteban et al. (1998) derived a factor of 
1.5 larger for ++ /H + from recombination lines than from CELs. Using several ions for 
a study of M8, Esteban et al. (1999) found again that the ionic abundances are always 
somewhat lower as derived from collisional forbidden lines than those from corresponding 
recombination lines. The mean abundance difference is a factor of ~2. Recently Tsamis 
et al. (2002) presented observations of 5 more H II regions - M 17 and NGC 3576 as well 
as the Magellanic Cloud H II regions 30 Doradus, LMC N11B and SMC N66. They found 
that the disparity was always in the same direction. For four of their objects, the ++ /H + 
abundance from O II recombination lines exceeded the corresponding value inferred from 
the nebular [O III] CELs by factors ranging from 1.8 - 2.7, while the factor was ~5 for 
LMC N11B. According to the limited statistics, apparently H II regions exhibit less of an 
abundance dichotomy than some of the PNs. 
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Most of the efforts to explain the abundance puzzle between collisional and recombina- 
tion values have attempted to do so by examining electron temperature (T e ) variations in 
the plasma. Because T e is not expected to vary dramatically within the (hydrogen) ionized 
region of a given nebula (e.g., Harrington et al. 1982; Baldwin et al. 1991; Rubin et al. 1991), 
a method developed by Peimbert (1967) has been extensively used. He developed a useful 
formalism by expressing the volume emissivity for a given spectral line in a Taylor series 
expansion around an average temperature defined such that the first order term vanishes. 
Since the [fractional] mean-square T e variation (called t 2 ) is expected to be small, only terms 
to second order need be retained. The resulting influence on elemental abundance deter- 
minations began with the work of Peimbert (1967), Rubin (1969), and Peimbert & Costero 
(1969). There is a vast literature now that measures the discrepancy between collisional and 
recombination abundances in terms of t 2 . 

The studies of H II regions mentioned above also derived the value of t 2 that forces the 
two abundance techniques to yield the same result. For the Orion Nebula, Esteban et al. 
(1998) required t 2 ~ 0.024 to reconcile the difference, while for M8, Esteban et al. (1999) 
found an even larger t 2 = 0.032 was needed on average for the various ions in common. 
Tsamis et al. (2002) found for 30 Doradus a similar value, t 2 ~ 0.03. For the drastically 
different ++ /H + abundances found for NGC 7009, Liu et al. (1995) discussed that it would 
be necessary to invoke t 2 ~ 0.1, which would then force agreement close to the higher 
recombination value - a value more than 2.5 times larger than the solar O/H of 7.41xl0~ 4 
(Grevesse & Sauval 1998). Such a large t 2 is not at all predicted by current theory/models 
(e.g., Kingdon & Ferland 1998). 

The current unsettled situation has led to efforts to broaden the study to include other 
variables besides T e to analyse the effects upon abundance determinations. One promis- 
ing avenue is to examine abundance derivation considering density variations, abundance 
variations, and T e variations in combination. Liu et al. (2000) took this approach in their 
investigation of the PN NGC 6153 with a two-phase empirical model. Pequignot et al. (2002) 
continued the study using photoionization models including two components with different 
heavy element abundances. 

This is an extension of an earlier paper (Rubin et al. 2002 - Paper I) that used HST 
STIS and WFPC2 observations to study the variation of T e in NGC 7009 as determined from 
the [O III] (4364/5008) flux ratio. Very low values for t 2 in the plane of the sky (t\) were 
found (always <0.01). In this paper, we focus solely on T e variations with the purpose to 
determine from the observational data the magnitude of t\ for the Orion Nebula. In section 
2, we present the HST observations and data reduction procedures. Section 3 includes a 
discussion and analysis of extinction. We determine the electron temperature distributions 
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in section 4. In section 5, we analyse the T e distributions in terms of average temperatures 
and fractional mean-square temperature variations in the plane of the sky. Section 6 provides 
a discussion and conclusions. 



2. HST OBSERVATIONS AND DATA REDUCTION 

The observations of NGC 1976 described here were taken as part of our HST Cycle 7 
program GO-7514. We observed with 4 different STIS long-slits: slits 1, 2, 4, and 5. These 
are shown in Figure 1. Our slit positions are chosen to cross several features, including an 
Herbig-Haro (HH) object, a proplyd, and the Orion Bar. Slit 1 passes through the position 
1SW which we observed with FOS and GHRS in Cycle 5 (see Figure 1 in Rubin et al. 1997). 
The centre of 1SW is at a, 5 = 5 h 35 m 14!71, -5°23'41'/5 (all positions are equinox J2000), 
18.5" S and 26.2" W of 6 l Ori C. This slit was also chosen to pass through proplyd P159- 
350 (O'Dell & Wen 1994). Slit 2 passes through position x2 observed with FOS in Cycle 
5 (see Figure 1 in Rubin et al. 1997). x2 is located in one of the most prominent arcs in 
the [O III] WFPC2 images. Slit 2 is parallel to slit 1, which allowed a substantial saving 
in overhead (OH). The position angle (PA) = 114.555°; the distance between slit centres is 
18.1". Note that the bottom of slit 1 and the top of slit 2 have a separation of only 0.211". 
The displacement between these two slits in the direction along the slit is 18.056". 

Slit 4 crosses the Orion Bar and is oriented to point toward 9 1 Ori C (the dominant 
exciting star), which also places it essentially orthogonal to the Bar. The southern tip passes 
through the Herbig-Haro object HH 203. Slit 5 passes through a very bright, sharply-defined 
"rim" of the Bar where a positional bifurcation begins. Slit 5 is parallel to slit 4 in order to 
take advantage of the same OH savings as is the case for slits 1 and 2. The PA = 139.068°; 
distance between centres of slits 4 and 5 is ~32". 

Observations were made as follows: data for slits 1 and 2 are from Visit 2 (1998 December 
7 UT), Visit 52 (2000 December 7) and Visit 72 (2001 December 16); data for slits 4 and 
5 are from Visit 5 (1998 December 22). All our data described here were acquired using 
the STIS/CCD with a 52" slit length and 0.5" slit width. Each visit comprises 2 orbits. 
Spectra were taken with gratings G430M (with wavelength settings: 3680, 4451, 4961) and 
G750M (settings: 5734, 6581). Each exposure was done in accumulation mode and at least 
two spectra were taken at each setting in order to cosmic-ray (CR) clean. The data sets 
we processed were those obtained after sufficient time elapsed from the observation dates 
in order that "best reference files" would be stable/finalized. We requested On-the-Fly 
Calibration for science files and Best Reference Files. After retrieving the data sets, we then 



- 5 - 



co-added and cosmic ray cleaned images using standard packages in IRAF. 8 Calibrations to 
produce 2-dimensional (2D) rectified images were then carried out. From these, we singled 
out specific emission lines for further investigation. Data for the [O III] 4364 A line as well as 
H7 4342 A were contained in the G430M/4451 grating setting; for the [O III] 5008, 4960 A 
lines and H/3 4863 A lines, G430M/4961 was used; for [N II] 5756 A, G750M/5734 was used; 
for Ha 6565 A and the [N II] 6550, 6585 A lines, G750M/6581 was used. All wavelengths in 
this paper are vacuum rest wavelengths. 

For G430M, the dispersion is 0.28 A per pixel for a point source and the plate scale 
0.05 arcsec/pixel (Leitherer et al. 2001, Chapter 13). For a uniformly filled slit with width 
0.5", a degradation in resolving power by a factor of 10 is expected to a spectral resolution 
of 2.8 A. For G750M, the dispersion is 0.56 A per pixel for a point source and the plate scale 
0.05 arcsec/pixel. For a uniformly filled slit with width 0.5", the spectral resolution would 
be 5.6 A. 

According to Paul Goudfrooij (private communication), the STIS CCD absolute flux 
calibration that is performed in the pipeline follows methods described in Instrument Science 
Report (ISR) 97-14 (Bohlin, Collins & Gonnella 1998) available through the STIS web site at 
www.stsci.edu/hst/stis/documents/isrs. That report considers only the L (low-resolution) 
modes specifically, but the flux calibration for the M (medium-resolution) modes is done 
the same way, i.e., by comparing the observed spectrum of a primary standard with a pure- 
hydrogen white dwarf model. The calibration observations were made with the 52" x 2" slit; 
thus relative transmission corrections are necessary to derive the absolute fluxes for the other 
slits (such as the 52" x 0.5" slit we used). This procedure is documented in ISR 98-20 (Bohlin 
& Hartig 1998). 

For the analysis presented here, we were interested mainly in the distribution of line 
flux along the slit spatial direction. This was accomplished with the IRAF routine blkavg in 
conjunction with specialized software tools that we developed ourselves. Even after applying 
the standard CR rejection there still remain many bad pixels due to CRs and/or hot pixels. 
There is considerable danger that including these can corrupt the flux values we seek. The 
program developed to eliminate these remaining bad pixels is called PIXHUNTER, which 
has been described briefly in an earlier paper (see Appendix A of Paper I). Once the columns 
containing the line have been cleaned for deviant pixels, we are ready to subtract an equiv- 
alent spectral range of continuum. We do this by using IRAF functions, including blkavg, 
to operate upon the appropriate sections of cleaned continuum. The 1-dimensional (ID) 
distribution of line flux versus spatial coordinate for the various emission lines of interest is 



8 IRAF is distributed by NOAO, which is operated by AURA, under cooperative agreement with NSF. 
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what we need for our subsequent physical analysis. 

We note that there is excellent agreement with a cross check of the ID results of flux 
versus spatial direction by comparing with ID results of flux versus wavelength for a cor- 
responding spatial sample. The latter were measured with the splot package. With this, 
the underlying continuum is fitted and the integrated line flux determined with the e-option 
(area under the line profile), which was preferable to fitting with a Gaussian profile. Because 
of the spectral impurity introduced by the relatively wide slit used, the line profiles have 
flatter tops and less extended bases (i.e., they are more "trapezoidal") than the Gaussian 
fits. It is also apparent that the Gaussian fit is overestimating the line flux. 

Both the [O III] 5008 and 4960 A lines were observed simultaneously with the 
G430M/4961 grating setting. Because both transitions arise from the same upper level, the 
intrinsic flux ratio depends only on the transition probabilities (A-values) and wavelengths. 
As reported previously for similar STIS observations of the planetary nebula NGC 7009 ob- 
tained under program GO-8114 (PI RR) (Paper I), what we found was a surprising variation 
in the F(5008)/F(4960) ratio with position along the slit. This amounts to a variation in 
the ratio of roughly 3.0±0.1. Furthermore, it appears more-or-less periodic with an ~3.5" 
cycle. As described in Paper I, according to Ted Gull (private communication), this is an 
instrumental effect and is a ratio of two fringe patterns. The source of the problem is a thin 
blocker filter that had to be matched with each grating and the best (and only) location that 
it could be placed was above the grating in a stable mounting. To attempt to do anything 
about fringing would probably require a dedicated HST/STIS calibration program. If there 
were fringing in the F(5008)/F(4364) ratio at the same level as for the F(5008)/F(4960) 
ratio, the ± 3.3% error would result in only a minor T e error, e.g., ± 100 K at T e = 10 4 K 
(see §4.1). 



3. EXTINCTION and REDDENING CORRECTION 

Before deriving the T e distribution from the STIS data, we first correct for extinction. 
This is calculated by comparing the observed F(Ho;)/F(H / 9) ratio with the theoretical ratio 
I(Ha)/I(H/3). We use a value of 2.88 assuming T e = 8500 K and N e = 5000 cm" 3 , Case B 
(Storey & Hummer 1995). 9 The extinction correction is done in terms of c(H/3), given by 



9 For conditions applicable to NGC 1976, we find using an online program (see Storey & Hummer 1995) 
that the I(Ha)/I(H/3) ratio used here will depart by less than 3% over the range 7500 < T e < 12500 K and 
10 3 <N e < 10 5 cm" 3 . 
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the relationship 

log[F(X)/F(H(3)\ = log[I(X)/I(H(3)\ - /(A) c(H(3), (1) 

where /(A) is the extinction curve. For the H7 (4342), 4364, 4863, 4960, 5008, 5756, 6550, 
6565, and 6585 A lines, the respective values for /(A) are 0.0856, 0.082, 0, -0.0153, -0.0225, 
-0.1237, -0.2185, -0.2202, and -0.2226, (Martin et al. 1996). This leads to, 

c(H(3) = 4.541 log[F{Ha)/F{H(3)} - 2.086. (2) 

The correction for extinction/reddening from observed to intrinsic flux for the 4364 and 
5008 lines is then given by, 

1(4364) = F(4364) 10 L082 c(H ^ ; 7(5008) = F(5008) 10 a9775 < m , (3) 

and for the 5756 and 6585 lines by, 

7(5756) = F(5756) 10 a8763 c{H(}) ; 7(6585) = F(6585) 10 a7774 c{m . (4) 

For the STIS data, we binned the pixels along the slit into tiles that are 0.5" square 
(matching the slit width). The STIS fiducial bars are excluded. This produced c(H/3) results 
for from 79 to 96 tiles, depending on the slit/visit observed, which are the same set used 
later for the T e analysis. The distributions of c(H/3) along the two slits (slit 1 and slit 2) 
observed in 3 separate visits (V2, V52, and V72) are remarkably similar. Furthermore, the 
distributions match well for the regions of slit 1 and slit 2 that are adjacent to each other. 
Recall that the separation between the bottom of slit 1 and the top of slit 2 is a mere 0.211" 
and that slit 2 is shifted 18.056" relative to slit 1 along the slit spatial direction. The average 
c(H/9) values without regard to any weighting for brightness are: for slit 1: 0.680, 0.699, and 
0.701 respectively for the three visits; for slit 2: 0.762, 0.798, and 0.766. For slit 4, it is 0.566 
and for slit 5, it is 0.611. We comment further on slit 1, which covers the proplyd P159-350 
and 1SW. Starting at the SE end of slit 1, the values are roughly flat at ~0.85. There is a 
spike to 1.43 in V2 at the tile containing most of the proplyd. The behavior is similar for 
V72 with a spike to 1.33 at the tile containing most of the proplyd. The adjacent nebula 
has for both visits a c(Rf3) value of ~0.9. Unfortunately, because P159-350 was centred in 
the East fiducial bar on V52, we cannot reliably measure c(R/3) there. Where slit 1 crosses 
1SW, c(H/9) is ~0.6, which agrees well with our previous spectroscopic value 0.605 at our 
position observed with the Faint Object Spectrograph (FOS-1SW, Rubin et al. 1998). At 
the end of slit 1 toward the NW, c(R(3) has decreased to ~0.4. 

We examined the vicinity of P 159-350 in more detail using the original resolution of 1 
pixel (0.05") for slit 1 in the spatial direction. The increase in flux and thus signal-to-noise 



- 8- 



(S/N) in the vicinity of the proplyd permits meaningful analysis here. For V2, we find that 
the peak observed flux occurs in pixel 311 for the Ha (5.64xl0 -11 erg cm -2 s _1 arcsec~ 2 ) 
R/3 (6.83xl(T 12 ), and 4364 (3.18xl0™ 13 ) lines, while the 5008 line reaches the highest peak 
(6.71xl0~ 12 ) at pixel 317 and a second relative maximum (5.14xl0~ 12 ) at pixel 305 with 
a clear trough between these, having a relative minimum (3.58xl0~ 12 ) occurring at pixel 
310. The F(Ha)/F(H/3) ratio peaks at 8.269 also at pixel 311 and hence the derived c(rL3) 
reaches a peak value there of 2.08. 

Repeating the analysis for V72, we find that the peak observed flux occurs in pixel 343 for 
the Ha (2.19xl0~ n ) and 4364 (8.11xl0~ 14 ) lines. The H/3 peaks in pixel 344 (3.19xl0~ 12 ) 
and is lower at pixel 343 (2.74x 10~ 12 ). The 5008 line reaches the highest peak (5.77xl0~ 12 ) 
at pixel 350 and a second relative maximum (4.58xl0~ 12 ) at pixel 336, again with a distinct 
trough between these, having a relative minimum (2.85 xl0~ 12 ) occurring at pixel 344. The 
F(Ha)/F(H / 9) ratio peaks at 7.994 at pixel 343 where the derived c(R(3) reaches a peak 
value of 2.01. The fact that the peak surface brightnesses are much higher for P159-350 in 
V2 compared with V72 is most likely because the brightest part of the proplyd was better 
sampled and/or aligned with the 1 pixel x 10 pixel (0.05" by 0.5") "smallest rectangular 
aperture" in the former visit than the latter. 



4. ELECTRON TEMPERATURE DETERMINATION 

4.1. [O III] Electron Temperature 

The electron temperature T e is derived from the intrinsic ratio I(5008)/I(4364) using 
the following relation, 

T e = 32966/[/n(/(5008)//(4364)) - 1.701] . (5) 

Effective collision strengths are from Burke, Lennon & Seaton (1989) for T e = 10 4 K. Transi- 
tion probabilities (A- values) are from Froese Fischer & Saha (1985). Note that this holds in 
the low-iV e limit, which should be valid for Orion where N e values are less than the critical 
densities (N crit ) for these lines. The lowest N crit ~6.4xl0 5 cm -3 for the 5008 line, which is 
well above N e values determined (e.g., Peimbert & Torres- Peimbert 1977; Osterbrock, Tran 
& Veilleux 1992; Esteban et al. 1998). 

We continue with the analysis using the tiles along the slit described above; the number 
of usable tiles varies with the slit/visit and ranges from 79 to 96. Equations (2), (3), and 
(5) are applied to the four emission lines to derive T e . Figure 2(a) shows the distribution of 
T e versus position along slit 1. For slit 1, position along the slit is relative to the location 
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of peak surface brightness (S.B.) in the Ha line, which occurs at P159-350. We set our 
zero corresponding to coordinates, a, 5 = 5 h 35 m 15!94, — 5°23'50"04 (Bob O'Dell, private 
communication). As mentioned in the last section, there is actually a relative minimum in 
the [O III] 5008 A S.B. close to that position. The lower curve in Figure 2(a) shows the 
observed 5008 S.B. (erg cm -2 s _1 arcsec -2 ), which is displayed (unsmoothed) at the pixel 
level providing 0.05" spatial resolution along the slit. The feature at the very left (SE) end 
of this slit is associated with the prominent arc, part of which we used to define our position 
x2 and which was better observed in slit 2. The open circles on the top (T e ) curve represent 
the individual tiles plotted at their midpoint. The dashed lines are a linear interpolation 
across the tiles that were deemed to have unreliable measurements because of proximity to 
the fiducial bars. There is a remarkably flat T e distribution with the notable exception of an 
upward spike to 13600 K at the tile containing the bulk of the proplyd Ha emission. For slit 1 
and V2, the behavior is similar with a sharp upward spike to 15050 K at the tile containing 
the bulk of the proplyd Ha emission. It is most unlikely that these high temperatures are 
realistic (these tiles are omitted in the statistical analysis to follow) and are the result of 
equation (5) not accounting for the higher electron densities, probably in excess of 10 6 cm -3 , 
in P159-350. Because the 5008 line will then suffer considerable collisional deexcitation, T e 
derived using equation (5) will be overestimated (e.g., Viegas & Clegg 1994). In the next 
section, we will evaluate the distribution of T e along slit 1, as well as along the other slits, 
in terms of T e variations. 

Figure 2(b) shows the distribution of T e versus position along slit 2, where position 
along the slit is relative to the location of the local peak S.B. in the 5008 line that occurs at 
x2. Using the position of P159-350 and the specified offset from slit 1 to slit 2 of +1.08s in 
RA and -8.15" in Dec, we find x2 coordinates: a, 5 = 5 h 35 m 16!96, -5°23'57"73. The lower 
relative maximum near offset +14" is where slit 2 passes through the "downstream" tail of 
P159-350. Again there is a remarkably flat T e distribution. 

Figure 3(a) plots the distribution of T e versus position along slit 4, where position along 
the slit is relative to the location of a feature in the bar seen in the 6585 S.B. distribution 
(Fig. 3(c)) discussed in the next section. Recall that this slit is aligned to point toward 
1 Ori C. While there is an underlying decrease in the 5008 S.B. with increasing distance 
from 9 1 Ori C, there is a leveling off (plateau) at the SE end of slit 4 of 5008 emission 
that remains substantial well beyond the bar. This is consistent with the WFPC2 images 
of the region; the composite "drizzled" image by Walsh (1998) that includes filter F502N 
(which covers the 5008 line) depicted in green, shows a green hue in a parabolic shape that 
appears to be a wake emanating from HH 203 and HH 204. Walsh's image is available 
at stecf.org/newsletter/webnewsl/orion/m42col_drizzle.jpg. The point was also made in an 
earlier discussion of HH 203 and HH 204 (O'Dell et al. 1997). Along this slit, the T e 
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distribution is also fairly flat, especially in the region of higher S/N (see §4.2). There is a 
clear increase in the amplitude of the T e fluctuations for the tiles SE of the bar due to poorer 
S/N in the 4364 A line. 

In Figure 3(b), we plot the distribution of T e versus position along slit 5. Here, the 
position along the slit is relative to the location of a feature seen in the 6585 S.B. distribution 
(Fig. 3(d)) discussed below. There is roughly a linear decrease in 5008 emission from the 
NW to the SE end of slit 5 with a "jump" occurring between 18-22" offset position. As was 
the case for slit 4, here too there remains substantial 5008 emission beyond the bar that is 
evident in the green hue at this location in the Walsh image. The T e distribution here is flat 
similar to that for slit 4; the increase in the T e amplitude range at the SE end of the slit is 
again due to poorer S/N in the 4364 A line. 

4.2. [N II] Electron Temperature 

The T e in the N + zone is derived from the intrinsic ratio 1(6585) /I(5756) (see equations 
(2) and (4)). Because the critical density for the 6585 line, ~7.7xl0 4 cm" 3 (at 10 4 K), is 
substantially less than N crit for the [O III] 5008 line (~6.4xl0 5 cm -3 ), we do consider various 
N e values when deriving T e . We also derive temperatures utilizing two different sets of N II 
effective collision strengths - those calculated by Lennon & Burke (1994) and by Stafford 
et al. (1994). We use the effective collision strengths for 10000 K; these do not vary much 
with the T e range of interest in Orion. The A-values used are discussed in Rubin et al. (1998, 
Appendix A with original references therein). 

For non-zero densities, our derivation of T e is done in an iterative fashion, starting 
with an initial estimate in the low-Ag limit. Then the volume emissivities (j values) for 
both the 5756 and 6585 lines are calculated solving the statistical equilibrium equations 
for the six lowest energy levels. The T e value is then recomputed using the intrinsic ratio 
1(6585) /I(5756), j(6585), and j(5756). The rapidly converging iteration is halted when T e 
changes by less than 1 K. 

In this paper, we are interested in assessing the amount of t 2 that can occur. Thus, 
we limit our analysis here by not including a detailed study of density variations, which is 
the subject of future papers by the present authors. Nevertheless, it is necessary to include 
some discussion of N e with regard to deriving the [N II] T e . As is well known, there is an 
inverse scaling of [N II] T e with N e . First we perform the calculations for T e in the low-Ae 
limit. This provides a firm upper limit to the T e distributions evaluated for the various 
slit/visit data. We also repeat the full computations for 4 other iV e 's: 1000, 2000, 5000, and 
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10000 cm~ 3 . The last value is more than high enough to bracket expectations for Orion, 
with the exception of objects like the proplyds as will be discussed below for P159-350. The 
value of 5000 is our best single N e to cover the region of slits 1 and 2 (e.g., Pogge, Owen & 
Atwood 1992). Their map of the electron density in the S + zone indicates that N e is lower 
in the regions of our slit 4 and slit 5; for these, we adopt a best single N e = 2000. It is these 
respective densities that are used to depict the distribution along slits 1, 2, 4, and 5 of [N II] 
T e in Figure 2(c), 2(d), 3(c), and 3(d) respectively. For the Figures, we have used the N + 
collision strengths of Lennon & Burke (1994) while results from this set, as well as Stafford 
et al. (1994) will be tabulated later. The lower curve in each of these panels shows the 6585 
surface brightness distribution along the slit. 

Figure 2(c) shows the tiled [N II] T e distribution versus position along slit 1 and is a 
positional match to Figure 2(a). The lower curve here shows the observed 6585 A S.B. at 
the unsmoothed pixel level. There again appears to be a relative minimum in the [N II] 6585 
S.B. at the centre pixel for P159-350 although it is not as well defined as the dip in the 5008 
S.B. there. As with the 5008 curve, the peak S.B. on the NW side of the proplyd exceeds 
that on the SE side. However, the overall 6585 emission is much narrower in P159-350 than 
is the 5008 emission. In totality, the structural behavior in the 6585 emission compared with 
that in 5008 is consistent with the N + region being more tightly confined than is the ++ 
region to the low- mass star of PI 59-350, which is expected if the proplyd is externally ionized 
by 9 1 Ori C. The prominent absolute maximum of the 6585 S.B. is very close to the 1SW 
position, 20.26" from P159-350, which is somewhat blocked by the NW-side fiducial bar. 
The narrow secondary peak at ~13" offset appears to be associated with the much broader 
secondary maximum near offset ~12" in the 5008 S.B. of Figure 2(a). There is perhaps an 
ionization boundary being observed here, although the slit orientation is far from ideal to test 
this. A slit pointing toward 6 1 Ori C would be better suited. Once again, there is a notably 
flat T e distribution with the exception of the proplyd. In Figure 2(c), we have omitted the 
tile with the highest derived T e of 20815 K. The highest T e in the plot is at 13000 K. Both 
these tiles are omitted from the subsequent statistical analysis. These high derived T e 's 
are undoubtedly pointing more to high N e values associated with the proplyd than high 
temperatures. The N e assumed for the calculation of T e in Figure 3(c) is 5000 cm -3 , while 
N e for P159-350 may exceed 10 6 cm -3 . Because the 6585 line emission will suffer enormously 
from collisional deexcitation at more realistic densities, the derived T e would be much lower. 
For slit 1, V2, there are two tiles at the proplyd that have an enormous upward spike in T e to 
29800 and 35850 K; both tiles are omitted from the statistical analysis in following sections. 

Figure 2(d) has the distribution of [N II] T e versus position along slit 2. There is a fairly 
flat T e distribution with some indication of a drop in T e SE of the lower fiducial bar. The 
maximum 6585 S.B. in the lower curve occurs near offset 28". This location is immediately 
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adjacent to the narrow peak at ~13" offset in Figure 2(c), just mentioned. The relatively 
small peak at x2 (offset 0) further supports the finding in the WFPC2 composite image that 
the emission arc which includes the position x2 is mainly an [O III] feature (see Figure 1). 
We note that x2 is part of the HH 529 complex and more specifically, the easternmost feature 
called 170-358 (see fig. 20 in Bally, O'Dell, & McCaughrean 2000). According to their study, 
including the kinematics, 170-358 is an expanding bow shock seen from the side. 

Figure 3(c) plots the distribution of [N II] T e versus position along slit 4. Position along 
this slit is relative to the location of the SE-most peak in S.B. in the 6585 line. We set our 
zero corresponding to coordinates: a, 5 = 5 h 35 m 21!34, — 5°24 / 48 / /21 determined from the 
position of our offset star 9 2 Ori A (HD37041) with a, 5 = 5 h 35 m 22!9, -5°24'57"9. There is 
a brighter feature that occurs where the slit crosses HH 203 at offset near —20". Comparison 
with Figure 3(a) shows no enhancement in 5008. This is consistent with the WFPC2 images 
of the region; the composite image in Figure 1 that includes filter F658N (isolating the 6585 
line) shows HH 203 as a reddish-orange object. There is a very flat T e distribution even over 
the wide range of S.B. and ionization structure along slit 4. 

In Figure 3(d), we plot the distribution of T e versus position along slit 5. Here, the zero 
position along the slit is the SE-most 6585 S.B. peak of where there is a bifurcation in the 
bar (see Figure 1). This feature appears to mark the boundary between the ionized region 
and the photodissociation region (PDR). Using the position of Q 2 Ori A and the specified 
offset from slit 4 to slit 5 of -1.87s in RA and -15.32" in Dec, we find: a, 5 = 5 h 35 m 19!74, 
— 5°25'08"14. Beyond this position, the T e distribution becomes less reliable because of 
poorer S/N in the 5756 A line. Overall, the T e is flat as has been the case throughout. 



5. FRACTIONAL MEAN-SQUARE TEMPERATURE VARIATIONS 

Our STIS analysis above presents results in the plane of the sky. The observations here 
do not address temperature fluctuation along the line of sight, which may be characterized in 
terms of the average temperature T and fractional mean-square T e variation (t 2 ) as defined 
by Peimbert (1967). 

/ T e N e N t dV 

T ° ~ J N e N z dV ' (6) 

/ (T e -T Q ) 2 N e N t dV 

T 2 jN e N t dV ' [<) 

where iVj is the ion density N(N + ) or N(0 ++ ). The integration in equations (6) and (7) 
is over the column defined by each tile, and along the line-of-sight (los). We are unable to 
measure the t 2 along the los for any column (cross section 1 tile). If there are t 2 along the 
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los, we can say that T(4364/5008) > T or T(5756/6585) > T (e.g., Peimbert 1967; Rubin 
et al. 1998). 

In the case of [0 III], for each tile, we have calculated T(4364/5008). Then the intrinsic 
flux 1(5008), fully correcting F(5008) for extinction (see equ. 3), in each tile is used in 
conjunction with T e = T(4364/5008) for that tile, and assumed constant along the los, to 
derive the following: 

/(5008) = 7^(5008) / N e NiT-°- 5 exp(-x/kT e )dl = K(5008,T e ) [ N e N i( il . (8) 



Here x is the excitation energy above the ground state for the upper level of the 5008 
transition, k is Boltzmann's constant, i^(5008) is known from atomic data, and i^(5008,T e ) 
has finally incorporated the known T e factor with the atomic constants. Here we again make 
the safe assumption of the low-N e limit (negligible collisional deexcitation) discussed earlier. 

For the case of [N II], for each tile, we have calculated T(5756/6585). Then the intrinsic 
flux 1(6585), fully correcting F(6585) for extinction (see equ. 4), in each tile is used in 
conjunction with T e = T(5756/6585) for that tile and the chosen value for N e , and assumed 
constant along the los, in the following relation: 

™ « S • < 9 > 

Here e(6585) is the normalized volume emissivity which is related to the usual volume emis- 
sivity j(6585) by e(6585) = j(6585)/(iV e iVj). e(6585) depends on the (fractional) population 
in a given level obtained by solving the 6-level atom for the specific T e and N e . It is not 
necessary to deal with the constant of proportionality for our purposes. 

Following Paper I, we define the average T e (T a) and fractional mean-square T e varia- 
tion (t\) in the plane of the sky. 

J J T e N e N t dldA 
T °' A ~ J J N e N t dldA ' (10) 

2 If (T e -T QA fN e N t dldA II T^NeNtdldA 

A " Tq,aI I N e NidldA Tl A I I N e N t dldA ' 1 } 

where dA represents an element of surface area in the plane of the sky and the integration 
over dl is for each tile along the los. The proper weighting for each tile is provided by either 
equation (8) or equation (9) depending on whether we are performing an [O III] or [N II] 
analysis. 
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The results using the STIS data for the various slits and visits to determine T 0:A and 
t\ are summarized in Table 1. These quantities for [O III] are entered in the first column. 
For slit 1, Visit 72 (see Figure 2(a)), T 0yA = 8258 K and t\ = 0.00925. We have excluded 
the spike in T e at the proplyd. The value for t\ is remarkably low, although not as low as 
the value of 0.0035 that we found for NGC 7009 (Paper I). The results from the other visits 
(V2 and V52) for slit 1 are close to those for V72. Our results for slit 2 and all three of 
the visits show T 0iA values in close accord and similar to those for slit 1. All these averages 
vary by less than 285 K. By comparing t\ for slits 1 and 2 in all 3 visits, we see that all 
six numbers are small, ranging from 0.00682 to 0.0129. For slits 4 and 5 the statistics are 
roughly comparable and differ from those for slits 1 and 2 in the following ways: T ^ A is 
several hundred K smaller and t\ is somewhat higher, although the largest value, 0.0176, is 
still notably small. A simple visual comparison of the T e plots in Figs. 2(a), 2(b), 3(a), and 
3(b) leads one to conclude that t\ looks larger for slits 4 and 5, particularly with the larger 
T e excursions in the regions SE of the bar. We note however that these data with poorer 
S/N have a minor effect on the statistics. There is natural biasing against these tiles due 
to weaker emission in both 4364 and 5008 as can be seen from equation (8), which in turn 
provides a smaller (J N e N { dl) -weighting of these tiles in equations (10) and (11). 

All the remaining columns in Table 1 pertain to the derivation of T ^ A and t\ for [N II]. 
These are displayed for 5 different electron densities and two different sets of effective collision 
strengths. As described, for slits 1 and 2, our preferred N e is 5000 cm -3 , while for slits 4 
and 5, it is 2000 cm -3 . We will now discuss the resulting numbers in these particular 
density columns and also using the left-side column with the Lennon & Burke (1994) collision 
strengths, unless stated otherwise. For slit 1, Visit 72 (see Figure 2(c)), T QtA = 10226 K and 
t\ = 0.00695. We exclude the high T e near P159-350 (see §4.2). The value for t\ is even 
lower than found for this slit /visit in the [O III] analysis; T QA is nearly 2000 K higher. There 
is close agreement among the N + T 0jA and t\ values for all the visits for slit 1; the above 
intercomparison between N + and ++ with the corresponding numbers in the respective 
visits is also true. 

For slit 2, the T ^ A numbers vary by more than 400 K between visits; the respective 
Tq ^'s are smaller than for slit 1 for V52 and V72 but slightly higher for V2. All of the t\ 
values are larger for slit 2 than for slit 1, reaching 0.0146. 

For slits 4 and 5, the N + Tq >a is several hundred to more than 1000 K smaller than for 
slits 1 and 2. Compared with the ++ T A , slits 4 and 5 are higher by 1343 and 1749 K, 
respectively. t A at 0.0175 is largest for slit 5 and nearly identical to the corresponding ++ 
value, which is also the highest in the first column. This larger t\ may be attributed to 
the larger excursions in T e in those tiles SE of the bar, where the S/N in the 5756 A line is 
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poorer. A comparison of Figure 3(c) and 3(d) clearly indicates a much lower 6585 S.B. SE of 
the bar for slit 5 than for slit 4. Again we emphasize that these data with poorer S/N play 
only a minor role in contributing to the integrals in equations (10) and (11). When there is 
weaker emission in both 5756 and 6585, as can be seen from equation (9), there is a smaller 
(J iV e Ni dZ)-weighting of these tiles in equations (10) and (11). 

The effect of substituting the Stafford et al. (1994) collision strengths is small for both 
slit 1 and slit 2 for all 3 visits: T 0j a decreases by less than 200 K; t\ decreases by roughly 
5-6 percent. The T ^ drop for slits 4 and 5 is slightly more, up to a 276 K difference; t 2 A 
decreases by ~6-7 percent. 

On the whole, the results presented in Table 1 allow us to reach two robust conclusions. 
First, t\ is always very small (< 0.0182) for either the ++ or N + analyses. Second, the T 0j a 
derived from the [O III] lines is always significantly less than T 0j a derived from the [N II] 
lines. These hold even when we consider for the N + analysis a wide range in N e that should 
amply bracket the great bulk of the plasma in Orion and also examine the influence of using 
alternative collision strengths. 

In Appendix A, we present an analysis of the uncertainties in the derived T e values as 
well as the t\ values. 

6. DISCUSSION AND CONCLUSIONS 

We extend our earlier study of electron temperature variations from the [O III] (4364/5008) 
flux ratio in the PN NGC 7009 (Paper I) to the Orion Nebula. In addition, we use our 
STIS long-slit observations to expand the analysis to measure T e variations from the [N II] 
(5756/6585) flux ratio. The observations here do not address T e fluctuation along the line 
of sight through the specific ++ region or, likewise, through the N + region. We analyse 
both the [O III] and [N II] data sets to derive the average T e and fractional mean-square 
T e variations in the plane-of-the-sky, which we call T 0>A and t\. We assume for each square 
column (projection of 1 STIS tile 0.5" square on the plane-of-the-sky) that the plasma along 
the line of sight is isothermal at the T(4364/5008) in the case of ++ or T(5756/6585) in 
the case of N + . For the latter case, we consider a large range of N e values, which should be 
sufficient for Orion, in order to produce Table 1. The analysis for each N e assumes that it 
is constant for the entire STIS slit length and throughout the sheet projected through the 
nebula along the los. 

Fluctuations in T e (and N e ) along the los are inevitable. We can make some com- 
ments about how our results for t A might be adjusted by T e variations along the los. The 
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relationship between T(4364/5008) and T for the ++ region is 



T(4364/5008) = T [1 + 0.5(91200/T C 
and between T(5756/6585) and T for the N + region is 



3)t 2 ] 



(12) 



T(5756/6585) = T [1 + 0.5(69000/T 



3)t 2 } 



(13) 



(e.g., Peimbert 1967; Rubin 1969). With the T(4364/5008)- or T(5756/6585)-values for 
Orion (see Figures 2 and 3) or indeed for H II regions in general, To will be smaller than 
these temperatures inferred from the forbidden line flux ratios here. Rather than repeat 
some examples of how the differences depend on t 2 , the reader is referred to our comments 
in the earlier paper (§7 in Paper I). 

We do not have the data here to characterize T e variations in 3-dimensions (3-D). It 
is useful to define an overall 3-D average T e (Toy) and fractional mean-square T e variation 
(ty). These single values apply for the entire source. Equations (6) and (7) define these 
specific values when the integration is over the entire volume. We note that for a spatially 
unresolved object (total integrated fluxes observed in the aperture), not the case for Orion, 
t v = t 2 and a calculation of t\ is meaningless. 

Our measurements of T e reported here are an average along each line of sight. Because 
each element of area treated in the plane of the sky represents a column which has already 
created a spatially averaged temperature along the los (e.g., see fig. 1 in Rubin 1969), it is 
likely that the value for t v is substantially higher than t\. Measurements of t 2 along various 
sight lines appear to be the most direct way to reliably gauge ty. Therefore, despite finding 
remarkably low t\, we cannot completely rule out much larger temperature fluctuations along 
the los. Further work, beyond the scope of this paper, is underway that will use modeling 
as well as additional observational data in an effort to better determine the relationship 
between t 2 , t\, and t v . 

Finding [N II] temperatures that are higher than [O III] temperatures is not a new result. 
For Orion, this has been known for years (e.g., Baldwin et al. 1991). What is new/significant 
here is that we have results from many more sight lines, more than 700 represented by the 
individual tiles, with improved spatial resolution (0.5 arcsec squares) and excellent spatial 
registration between the [N II] and [O III] data sets that HST affords. Of the roughly 700 
tiles, about half represent independent sight lines, with the other half having spatial overlap 
due to the repeated visits for slits 1 and 2. 



The fact that the [N II] temperatures are higher is also expected on theoretical grounds 
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and again, not something novel that we have uncovered. We enumerate three factors that 
contribute to the predicted inequality. First, the cooling of the [O III] 5008, 4960 A lines in 
the ++ region is uniquely efficient for nebular conditions and elemental abundances that 
prevail in Galactic H II regions, including Orion. On the other hand, in the "singly ionized" 
region (where the dominant O and N ions are + , N + ) there is no coolant that is nearly as 
efficient as [O III] 5008, 4960 is in the ++ zone. Additionally for Orion, the blister geometry 
likely results in higher average electron densities in the singly ionized region compared with 
the ++ region because the former is closer to the PDR. The higher N e in the N + region 
would also contribute to less efficient cooling. 

Second, there is the predicted hardening of the stellar ionizing photons at progressively 
larger distances from the exciting star. This is due to the functional form of the H photoion- 
ization (predominantly) cross section, which diminishes steeply with higher frequency from 
its value near threshold. At larger distances from the exciting source, the average energy 
per photoionization will increase and thus the heating rate will also increase. This causes a 
rise in T e with an increase in distance (other factors being equal). This has been known for 
many years (e.g., Rubin 1968). 

Third, there is a possibility that there is a contribution to the production of [N II] 5756 A 
emission (and to a lesser extent 6585) by recombination and cascading. Under some condi- 
tions, this may provide significant routes into the upper energy level of the 5756 transition 
that are not negligible compared with collisional processes. This was examined by Rubin 
(1986). If there were a significant "recombination" contribution to the observed emission, 
we would have overestimated the T e derived from the [N II] 5756/6585 flux ratio. With our 
two independent, detailed models for the Orion Nebula (Baldwin et al. 1991 and Rubin et al. 
1991), the recombination contribution appears to be small. We discussed this in a subse- 
quent paper where these models were "retrofitted" with more in-common input parameters, 
including the Stafford et al. (1994) N + collision strengths, to allow closer intercomparisons 
(see end of §4 in Rubin et al. 1998). 

All three of the above effects are included in our photoionization, plasma simulation 
models just mentioned. The predicted [N II] and [O III] temperatures for the entire volume 
are 8721 K and 7704 K in the Rubin et al. (1991) model and 8649 K and 7692 K in the retrofit 
(Rubin et al. 1998) using Stafford et al. cross sections. When the Lennon & Burke (1994) 
N + collision strengths are used instead of the Stafford et al. set, the respective temperatures 
become 8706 K and 7701 K. 

Finally, we encapsulate the findings here for the proplyd P159-350. We find large local 
extinction as evidenced by the dramatic increase in the observed F(Ha)/F(H/3) ratio along 
slit 1. For V2, this ratio peaks at 8.27 which implies a c(Hj3) of 2.08 from the Balmer 
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decrement. Similar values are found for our V72 observations: F(Ha)/F(H/3) peaks sharply 
at the proplyd position reaching a ratio of 7.99, where the derived c(H/3) = 2.01. Because 
the adjacent nebular values for c(H/3) are much lower, the extinction must be associated 
with the P 159-350 environs. 

A comparison of the [O III] 5008 and [N II] 6585 surface brightnesses in the vicinity of 
P159-350 (see Figures 2(a) and 2(c)) shows for both lines that the S.B. on the NW side of 
the proplyd exceeds that on the SE side. However, the 6585 emission is much narrower than 
is the 5008 emission which provides evidence of ionization stratification consistent with the 
N + region being more tightly confined than is the ++ region to the low-mass star of PI 59- 
350. This "inverse H II region" behavior is precisely what is expected due to the external 
ionization source 6 1 Ori C. 

The derived T e distributions were all notably flat except at the position of the proplyd 
observed in slit 1. The very high derived T e 's at the location of P159-350 are no doubt due 
more to high N e values associated with the proplyd than high temperatures. As stated, the 
calculations used here to assess T e have not accounted for N e values likely to exceed 10 6 cm~ 3 
in P159-350. Because both the 5008 line and even more so the 6585 line emission will suffer 
substantial collisional deexcitation at such high densities, the derived T e for both [O III] and 
[N II] would be much lower. 

Each of the above mentioned facets of P159-350 deserves further study. The smallest 
spatial resolution element was 1 pixel (0.05") in the spatial direction by 10 pixels (0.5") in 
the dispersion direction. It would be highly desirable to obtain higher spatial resolution 
of PI 59-350, which can be readily achieved by observing with a narrower STIS slit. The 
narrower slit would enable better spectral resolution as well and permit a detailed mapping 
of P159-350 in the C III] 1907, 1909 A lines. The ratio of their fluxes is an excellent diagnostic 
of N e even at the high N e values expected in proplyds. Improved spectral resolution in the 
optical lines will also help sort out kinematics from spatial structure, which was a difficulty 
with the 0.5" slit width we used. 

This paper is based on observations made with the NASA/ESA Hubble Space Telescope, 
obtained at the Space Telescope Science Institute (ST Scl). Support for proposal #GO- 
7514 was provided by NASA through a grant from ST Scl. ST Scl is operated by the 
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Bob O'Dell, Chris Ortiz, Aaron Svoboda, and Abby Wong. We thank Bob O'Dell and Janet 
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A. Error Analysis 

The values found for t\ are very small. It is important to determine whether even 
these small fluctuations in T e from tile to tile are a real signal or the result of measurement 
uncertainty. 

From equation (5) for ++ , it can be shown that the mean square uncertainty in tem- 
perature, expressed as t^, has the following dependence on the uncertainties 8I n and 8I a of 
the nebular (I n ) and auroral (I a ) line fluxes, respectively: 



t 2 e = (ST e /T e ) 2 = % 



5L 



(Al) 



where C\ = 32966. We use this as an approximation for N + too, with c\ = 24933. 



The uncertainty 51 is principally the measurement error in the line flux for a tile. The 
S/N for the flux of the auroral lines ([O III] 4364 and [N II] 5756 here) is much lower than 
that for the nebular lines ([O III] 5008 and [N II] 6585 here), and so the uncertainty of the 
weaker auroral line will dominate the uncertainty 5T e . The uncertainty 51 also depends on 
the (differential) extinction correction, but in the present case this is a small effect compared 
to the uncertainty in the observed auroral line fluxes (hence, I/5I ~ F/5F = S/N ). 

The three independent methods used to calculate the S/N of the auroral lines are dis- 
cussed in the following sections. We use T e ~ T 0)J 4 to calculate the corresponding t^. 



A.l. Method 1 

The flux integrated over the line profile as found using the "blkavg sum" technique in 
IRAF is 



m 



^ ^ p pline+cont ^ ^ p^pcont (A2) 

i=l U j=l 

where Fi is the monochromatic flux for the tile (erg cm -2 s _1 A -1 ), D is the dispersion 
conversion (A/pixel), n is the number of pixels in the spectrum used to define the continuum 
and m is the number of pixels used for the line + cont. Since the auroral line is weak, the rms 
fluctuation in the line + cont can be assumed to be the same as in the continuum, 5Fj ont . 
Thus, the rms of the line flux is 

5 piine = D^{5F™ nt )^l + m~/n~. (A3) 
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In our measurements we used m = n. The above three equations are used to create the 
entries for method 1 (Ml) in Table 2. Values are given for both the [O III] 4364 and [N II] 
5756 lines for all slits (for slits with multiple visits, a representative case is given). 

A. 2. Method 2 

In this method, the line profiles from the STIS spectra are fit with a template (a Gaussian 
convolved with a slit of width 0.5") using a combined Gauss-Newton and modified Newton 
algorithm (Numerical Algorithms Group (NAG) - E04FDF) to find the minimum least- 
squares solution and the corresponding variances of the variables used in the fit (NAG- 
E04YCF). The slit width is held constant, but due to anamorphic magnification in the 
dispersion direction, the plate scale differs from line to line which results in different slit 
widths expressed as pixels (see Bowers & Baum 1998). 

The variables used in the template fitting are the flux F line , the FWHM and central 
wavelength of the Gaussian, and the slope and mean brightness of the continuum. The 
variances reflect how well the model fits the data. Keeping all of these as variables gives one 
a conservative estimate of the variance of each. Note that a poor template would lead to an 
overestimate of the error. The templates were tested on the strong nebular lines and generally 
fit very well. However, there are instances where the surface brightness changes across the 
slit (in the dispersion direction), so that the implicit assumption of uniform illumination is 
not perfect. Therefore, the estimated uncertainty from method 2 is an upper limit. As is 
seen in Table 2 it agrees well with the error estimated using method 1. 

A.3. Method 3 

This method utilizes the multiple sets of observations for slit 1 and slit 2 in the three 
visits V2, V52 and V72 and examines the reproducibility of the line fluxes and T e . There 
are adequate fiducial sources visible in all three visits so that the spectra may be shifted 
in the spatial direction to achieve alignment to ~1 pixel. The proplyd P159-350 in slit 1 is 
particularly useful for this and serves for slit 2 as well because the relative positional HST 
offsets between slit 1 and slit 2 are not subject to the acquisition uncertainty. There is 
also positional alignment uncertainty in the dispersion direction which is more difficult to 
estimate. For all three visits, P159-350 was seen in the 0.5" wide slit; although it was blocked 
for the most part by the East fiducial bar in V52, some of the emission still peeks out both 
ends. We noted at the end of section 3 that P 159-350 has a higher peak surface brightness 
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in V2 than in V72. We estimate that the positional alignment of the tiles perpendicular 
to the spatial direction between the three visits is within 0.2" so that there is substantial 
area overlap (perhaps 60%) of our tiles between the three visits. Comparison of the nebular 
line strengths from visit to visit shows excellent reproducibility (< 3% per tile for the ++ 
line). This is much smaller than the auroral line differences found among successive visits. 
The measurement error of the weaker auroral lines is of course higher. Furthermore lack of 
registration is a more serious issue for the auroral lines (from a much higher energy level) 
if there are real T e fluctuations on this spatial scale; in the limit of no overlap, one would 
expect t 2 , ~ t\ even if there were no measurement errors. Therefore, having t\ from method 3 
larger than from Ml and M2 is not unexpected. 

For each tile in common between V2, V52 and V72, both an average value (F a ) and 
the standard deviation (5F a ) of any single determination of F a were calculated. The repre- 
sentative value entered as M3a in Table 2 is the median SF a /F a for each of the slits; as in 
methods 1 and 2, the 5T e listed is what is implied by this S/N, through equation Al. 

For each in-common tile, we also calculated T e for each visit, then the average T e and the 
actual standard deviation (5T e ) of any single determination of T e . The median ST e for each 
slit is included in Table 2 as M3b. From these median values we calculated the tabulated 
t\. Methods 3a and 3b agree closely, as they should given the good reproducibility of the 
nebular lines (i.e., the effect of 5F n in equation Al is negligible). 

The analysis summarized in Table 2 shows that the measured t\ is a real signal not 
dominated by the noise and/or measurement errors. The two exceptions are for the ++ 
measurements for slits 4 and 5 which include tiles beyond the Orion bar (SE part of the slit) 
where the surface brightness for lines of that ion, and hence the S/N, is low. Accordingly, 
we repeated the entire analysis for only the NW half of these slits. The results recorded in 
Table 2 show a lower t A , as anticipated, and a lower t 2 e as well which is a smaller fraction of 
i A . 

Where there is a contribution of measurement error to the apparent t A , the actual t\ is 
smaller, by about t\. This reinforces our conclusion that t\ is very small. 
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Table I. Values of T 0) a and t\ 



o++ 



N+ 



N e (cm- 3 ) 



1000 



2000 



5000 



10000 



Slit 1, Visit 2 
8151 b 
1.04 b 

Slit 1, Visit 52 
8232 
0.887 

Slit 1, Visit 72 
8258 
0.925 

Slit 2, Visit 2 
8142 
1.15 

Slit 2, Visit 52 
8074 
1.29 

Slit 2, Visit 72 
8358 
0.682 

Slit 4, Visit 5 
7790 
1.57 

Slit 5, Visit 5 
7678 
1.76 



11104 
0.860 



11236 
0.647 



11232 
0.765 



11171 
1.05 



10675 
1.60 



10909 
1.43 



9473 
0.985 



9789 
1.82 



(10605) c 
(0.795) 



(10727) 
(0.592) 



(10723) 
(0.703) 

(10665) 
(0.956) 

(10208) 
(1.48) 

(10424) 
(1.31) 

(9107) 
(0.918) 

(9394) 
(1.69) 



10881 
0.844 



11010 
0.633 



11005 
0.750 



10945 
1.02 



10463 
1.57 



10691 
1.40 



9298 
0.968 



9603 
1.79 



(10458) 
(0.785) 

(10578) 
(0.584) 



(10573) 
(0.694) 

(10517) 
(0.942) 

(10068) 
(1.46) 

(10280) 
(1.30) 

(8991) 
(0.906) 



(9271) 
(1.66) 



10670 (10317) 
0.831 (0.777) 



10114 (9929) 
0.791 (0.752) 



9381 (9386) 
0.730 (0.712) 



10796 (10435) 10232 (10042) 9488 (9491) 

0.620 (0.576) 0.584 (0.553) 0.531 (0.517) 



10791 (10429) 10226 (10036) 9481 (9485) 

0.736 (0.685) 0.695 (0.659) 0.636 (0.620) 

10732 (10373) 10170 (9982) 9429 (9433) 

1.00 (0.929) 0.941 (0.890) 0.853 (0.831) 

10262 (9933) 9733 (9564) 9038 (9046) 

1.54 (1.44) 1.46 (1.39) 1.33 (1.30) 

10484 (10141) 9938 (9760) 9221 (9227) 

1.37 (1.28) 1.30 (1.23) 1.18 (1.15) 

9133 (8879) 8698 (8572) 8125 (8143) 

0.952 (0.896) 0.906 (0.867) 0.837 (0.821) 



9427 (9151) 8962 (8825) 8353 (8369) 

1.75 (1.64) 1.66 (1.58) 1.52 (1.49) 



a clectron density assumed for N+ emission-line region; denotes the low density limit 
b upper row T 0j a (K); lower 100 t\ 

c using effective collision strengths from Lcnnon & Burke 1994 (Stafford ct al. 1994) 
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Table 2. Weak line uncertainties and the associated fractional mean-square temperature 

variation, t\ 



T ilTIP 


Slit 


Visit 

V lol u 


MpthnH a 


(SF /F ) j- 


T 

1 0,A 


ST 


1 00/ 2 


1 00/ 2 , 


4364 


1 


V72 


Ml 


0.1434 


8258 


297 


0.129 


0.925 








M2 


0.1132 




234 


0.080 








b 


M3a 


0.1826 




378 


0.209 








b 


M3b 






388 


0.221 




4364 


2 


V72 


Ml 


0.2003 


8358 


424 


0.258 


0.682 








M2 


0.2577 




546 


0.427 








b 


M3a 


0.2593 




550 


0.432 








b 


M3b 






513 


0.377 




4364 


4 


V5 


Ml 


0.3129 


7790 


576 


0.547 


1.57 








M2 


0.5350 




985 


1.60 




4364 


4NW 


V5 


Ml 


0.2435 


7665 


434 


0.320 


0.973 








M2 


0.1550 




276 


0.130 




4364 


5 


V5 


Ml 


0.4274 


7678 


764 


0.991 


1.76 








M2 


0.4527 




810 


1.11 




4364 


5NW 


V5 


Ml 


0.3235 


7683 


579 


0.568 


1.42 








M2 


0.3310 




593 


0.595 




5756 


1 


V72 


Ml 


0.0499 


10226 c 


209 


0.042 


0.695 c 








M2 


0.0616 




258 


0.064 








b 


M3a 


0.0937 




393 


0.148 








b 


M3b 






388 


0.144 




5756 


2 


V72 


Ml 


0.1188 


9938 c 


471 


0.224 


1.30 c 








M2 


0.1090 




432 


0.189 








b 


M3a 


0.1879 




744 


0.561 








b 


M3b 






729 


0.538 




5756 


4 


V5 


Ml 


0.1491 


9133 d 


499 


0.298 


0.952 d 








M2 


0.1226 




410 


0.202 




5756 


5 


V5 


Ml 


0.1459 


9427 d 


520 


0.304 


1.75 d 








M2 


0.1284 




458 


0.236 





a Methods 1-3 (M1-M3) are discussed in the text. 
b V2,V52 and V72 are used for Method 3. 

c N e = 5000 cm' 3 using effective collision strengths from Lennon & Burke 1994 
d N e = 2000 cm" 3 using effective collision strengths from Lennon & Burke 1994 
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Fig. 1. — This shows the position of the four STIS slits overlaid on a composite WFPC2 
image with [O III] (5008 A) blue, Ha green, and [N II] (6585 A) red. The CCD observations 
presented here are taken with the STIS long-slit 52" x 0.5". The middle portion of each slit 
shown in black represents the smaller length (28") for which there are other data obtained 
with the MAMA detectors. Abbreviated labels stand for positions described in the text: P 
for the proplyd P159-350; H3 for HH 203. Note that the star 6 2 Ori A is the unlabeled red 
feature just above the 4 in Slit 4. N is up and E to the left. 
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Fig. 2.— (a) Slit 1, V72: Plot of T e determined from the [O III] 4364/5008 flux ratio ver- 
sus position along the STIS long-slit. The analysis is in terms of tiles that are 0.5" square 
(matching the slit width). The open circles represent the individual tiles plotted at their 
midpoint. The dashed straight lines are interpolations across the two fiducial bars, indicated 
by the gray area, where data are unreliable. Positional measurement along the slit is from 
~SE to ~NW (see Figure 1). The zero point is described in the text. The bottom curve 
shows the observed [O III] 5008 A surface brightness in units of erg cm -2 s _1 arcsec -2 
displayed unsmoothed at the pixel level. (b) same as (a) for Slit 2, V72. 

(c) same as (a) except the upper curve is a plot of T e determined from the [N II] 5756/6585 
flux ratio assuming N e = 5000 cm -3 ; the lower curve shows the observed [N II] 6585 A 
surface brightness. (d) same as (c) for Slit 2, V72. 
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Fig. 3. — (a) Slit 4, V5: Same as Fig. 2 (a) (b) same as (a) for Slit 5, V5. 

(c) Slit 4, V5: Same as Fig.2 (c) except that N e = 2000 cm -3 is assumed. 

(d) same as (c) for Slit 5, V5. 



